On the RKKY range function of a one dimensional non interacting electron gas 
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We show that the pitfalls encountered in earlier calculations of the RKKY range function for a 
non interacting one dimensional electron gas at zero temperature can be unraveled and successfully 
dealt with through a proper handling of the impurity potential. 
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The apparently straightforward evaluation of the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) range func- 
tion, or more generally, of the linear density modula- 
tion 5n(x) induced at zero temperature in a non inter- 
acting one dimensional electron gas by a localized static 
impurity modeled with a delta function potential has 
proven surprisingly troublesome. While the original cal- 
culation gave an incorrect answer more recent investi- 
gations appear to suggest that only a certain procedure, 
known to lead to a physically sensible answer, should be 
employed |2|- 

The original and most popular procedure is based on 
the standard theory of linear response Q. If the impurity 
potential is assumed to be of the form U(x) = ^^■S(x) 
(where u is a suitable wave vector) then one has 
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where L is the system length and Xo(<7>0) is the static 
Lindhard response function in one dimension given by 
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where the notation V stipulates that the principal part of 
the sum must be taken. Substituting (J2J) in JQ) and using 
the zero temperature occupation numbers immediately 
leads to the formula 
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The interpretation and the handling of this expression 
rest at the origin of the problem at hand. 

If one simply proceeds to explicitly evaluate the inte- 
gral over k in Eq. ©, or, which is the same, that of 
the result is well known and is given by 
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At this point the second integral over q can be readily 
evaluated leading to 



5n{x) 
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where the sine integral function appears This result 
provides the correct answer to the problem. In particular 
the large distance behavior turns out to be 
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an expression displaying the expected decay and Friedel 
oscillations p|- One should notice at this point that, in 
view of the singular behavior of the integrand at the ori- 
gin, however tempting, the order of the k and q integra- 
tions in Eq. J2J cannot be freely exchanged 2J . Doing so 
leads to the manifestly unphysical answer[l| 
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This pitfall is unique to the one dimensional case for in 
two and three dimensions formulas equivalent to © can 
be derived and safely evaluated by exchanging at will the 
order of the intervening integrations. 

An alternative, appealing and equally physically valid 
procedure to obtain Sn(x) is in the case offered by ex- 
tracting the leading linear term of the formula 

n(x)-n = 2 J2 (\Mx)\ 2 -\M*)\ 2 ) > (8) 

\k\<k F 

valid for single Slater determinant states. Now, if the 
usual first order perturbation theory result for 4>k(x) is 
used in ©, one is immediately led to an equation that 
differs from merely for the exchange of the order of 
the quadratures lj. Since, as observed above, such an 
expression does lead to an unphysical result, it has been 
suggested that this route is unphysical and should be 
avoided 2j. We find that this conclusion is unwarranted 
for, as we show next, the difficulty lies here with the use 
of perturbation theory, which is invalid, and not with 
Eq. iJSJl per se. 

To prove our assertion we observe that the exact delo- 
calized eigenstates of the Schrodinger equation 
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can be written as a superposition of the two following leading to 
(normalized) scattering states 
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where fc is limited to positive values only. By taking the 
modulus square of the sum of ip^ (x) and ipk- (x) and 
by summing over k > in © one obtains the following 
exact expression for the electronic density n(x): 
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At this point we are left with extracting the linear term 
in u. This must be done with some care. In particular 
one must resist the temptation to simply drop the last 
term in Eq. (|12l) and, at the same time, to neglect the u 2 
in the denominator of the second one. To do so coincides 
with making use of first order perturbation theory and 
leads exactly to the original unphysical result of Eq. J7J . 
While it is safe to handle the second term of i|12[l as just 
described, the third term does contribute a first order 
term that can be readily extracted by making use of the 
relation 
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a result that can be readily seen to coincide with the 
correct answer of Eq. JSJ . 

We have therefore shown that a proper handling of the 
impurity potential allows one to correctly carry out the 
calculation of the linear density modulation via either the 
response function method of Eq. or the alternative 
direct procedure offered by ©P- 

We conclude by remarking that a proper treatment of 
the problem of the effects of a static localized impurity 
in an interacting one dimensional electron liquid can be 
achieved by means of the Luttinger liquid model. The 
problem is highly non trivial as, while Friedel-like oscil- 
lations with an amplitude decay generally ruled by the 
interaction strength exist at intermediate distances [3j, for 
large distances the physics of the phenomenon is non lin- 
ear in the impurity potential 7] . Finally a recent dis- 
cussion of the effects of the Coulomb interaction on the 
Friedel oscillations of an electron liquid in two and three 
dimensions can be found in reference 0. 

Acknowledgments 



The authors wish to thank A. W. Overhauser and 
George E. Simion for useful conversations. 



[1] C. Kittel, in Solid State Physics, vol. 22 (Academic Press, 

New York, 1968), edited by F. Seitz, D. Turnbull, and H. 

Ehreinreich (see also erratum ibidem). 
[2] Y. Yafet, Phys. Rev. B 36, 3948 (1987). 
[3] G. F. Giuliani and G. Vignale, Quantum Theory of the 

Electron Liquid (Cambridge University Press, Cambridge, 

2005). 

[4] I. Gradshteyn and I. Ryzhik, Table of Integrals, Series and 

Products (Academic Press, San Francisco, 1965). 
[5] J. Friedel, Advanc. Phys. 3, 446 (1954). 



[6] Of course the correct answer can also be arrived at by 
calculating 8n(x) at finite temperature and then taking 
the T — > as done for instance in V. Litvinov and V. K. 
Dugaev, Phys. Rev. B 58, 3584 (1998). This sheds no light 
onto the problem at hand. 

[7] R. Egger and H. Grabert, Phys. Rev. Letters 75, 3505 
(1995). 

[8] G. Simion and G. F. Giuliani, to be published. 



